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ABSTRACT 

The broad-line radio galaxy 3C120 is a powerful source of both X-ray and radio 
emission including superluminal jet outflows. We report on our reanalysis of 160 ks of 
Suzaku data taken in 2006, previously examined by Kataoka et al. (2007). Spectral fits 
to the XIS and HXD/PIN data over a range of 0.7-45 keV reveal a well-defined iron K 
line complex with a narrow Ka core and relativistically broadened features consistent 
with emission from the inner regions of the accretion disk. Furthermore, the inner 
region of the disk appears to be truncated with an inner radius of rj n = 11.7jJ' 2 r g . If 
we assume that fluorescent iron line features terminate at the inner-most stable circular 
orbit (ISCO), we measure a black hole spin of a < —0.1 at a 90% confidence level. A 
rapidly spinning prograde black hole (a > 0.8) can be ruled out at the 99% confidence 
level. Alternatively, the disk may be truncated well outside of the ISCO of a rapid 
prograde hole. The most compelling scenario is the possibility that the inner regions of 
the disk were destroyed/ejected by catastrophic instabilities just prior to the time these 
observations were made. 

Subject headings: galaxies: active — galaxies: individual (3C120) — galaxies: nuclei - 
galaxies: Seyfert — X-rays: galaxies 



1. Introduction 



The central engines of active galactic nuclei (AGN) are powered by the accretion of material 
onto a supermassive black hole (SMBH). Some fraction of the AGN population, the radio-loud AGN, 
produce powerful, relativistic jets (Donoso, Best, & Kauffmann 2009) — but the physical processes 
that distinguish radio- loud from radio-quiet AGN remain mysterious. It has been suggested that 
the spin of the SMBH plays an important role in determining whether an AGN is radio-loud 
or radio-quiet, a hypothesis known as the "Spin Paradigm" (Sikora et al. 2007). It is commonly 
accepted that relativistic jets are powered by the extraction of rotational energy from the ergoregion 
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of a spinning black hole (Blandford & Znajek 1977; McKinney, Tchekhovskoy, & Blandford 2012). 
However, whether the spin determines the radio-loudness of an AGN (i.e. whether there is an 
additional "clutch" that switches on/off the transmission of spin power into the jet) is still very 
much up for debate. 

The spin of a black hole with angular momentum J is described by the dimensionless spin- 
parameter a = cJ/GM 2 where \a\ < 1. In order to test the spin paradigm it is important that we 
measure the spin parameter in both radio-quiet and radio- loud sources. Previous measurements of 
several radio-quiet objects have found the central SMBH to be in a state of rapid prograde spin 
(Brenneman et al. 2011); this already presents a challenge to the spin paradigm. Our goal is to 
address the other side of this problem by studying the central engine structure and black hole spin 
of a radio- loud object. 

The current best technique for measuring SMBH spin is by studying the relativistically broad- 
ened spectral features produced when the surface of the accretion disk is irradiated by hard X-rays. 
The most prominent of these reflection spectrum features is the fluorescent Fe-Kck emission line 
found at 6.4 keV in the rest frame. The line emanates from the inner most regions of the accretion 
disk where it is broadened and skewed towards low energies by gravitational redshift and relativistic 
Doppler shift (Fabian et al. 1989). The line is clearly visible against the background continuum 
spectrum and serves as the strongest diagnostic for gravitational physics around SMBHs currently 
known (see reviews by Reynolds & Nowak 2003; Miller 2007). If one assumes that the line emitting 
region of an accretion disk extends down to the inner-most stable circular orbit (ISCO) (Reynolds 
& Fabian 2008) then, by modeling the profile of the broad FE-Ka line, the location of the ISCO can 
be discerned in units of r g = GM/c 2 , thereby giving us the spin (Brenneman & Reynolds, 2006). 

The broad-line radio galaxy (BLRG) 3C120 (z=0.033) is a prime candidate for such a study. 
It is the brightest BLRG in the X-ray sky and possesses a powerful one-sided jet on 100 kpc scales 
(Walker et al. 1987, Harris et al. 2004) with observed superluminal motion (f3 a pp = 8.1) indicating 
a jet inclination of i < 14° (Zensus 1989, Eracleous & Halpern 1998). Early X-ray observations by 
ASCA detected an extremely broad iron Ka line (Grandi et al. 1997, Reynolds 1997, Sambruna 
et al. 1999). A significantly narrower line complex was found in later observations by RXTE 
(Eracleous, et al. 2000, Gliozzi et al. 2003) and BeppoSAX (Zdziarski k Grandi 2001). Analyses 
of XMM-Newton data produced mixed results, with Ballantyne, Fabian & Iwasawa (2004) finding a 
weak broad iron line (albeit with a poorly constrained inner emission radius) but Ogle et al. (2005) 
claiming no relativistic line signatures. Aided by the hard X-ray coverage, observations of 3C120 
by Suzaku found both narrow iron emission lines and features consistent with a relativistically 
broadened line (Kataoka et al. 2007; hereafter K07). K07 were able to constrain the inner edge of 
the line emitting region to be r- m = 8.6^Q gr 9 . While not commented upon in K07, this begs the 
question as to whether 3C120 may be driven by a rapidly retrograde spinning SMBH (whose ISCO 
is at 9r g ). 

In this Letter, we present a reanalysis of the 2006 Suzaku data on 3C120 using updated in- 
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strument calibrations and the latest models/methodology for studying relativistically broadened 
reflection features. Section 2 contains an outline of the observations and data reduction steps 
taken. Section 3 presents the spectral analysis which fundamentally reconfirms the results of K07, 
and Section 4 discusses the implications of our analysis. Our results indicate either a retrograde 
spinning black hole, or a disk that is truncated beyond its ISCO. We suggest a possible relation to 
the disk-jet cycles observed to occur in BLRGs such as 3C120. 

2. Observations and Data Reduction 

The radio galaxy 3C120 was observed with Suzaku starting on 2006 February 7 beginning at 
03:20 (UT) and ending on 2006 March 03 20:29 (UT). A total of four quasi-continuous observations 
were made, each with approximately 40 ks of good on-source exposure. Suzaku is equipped with four 
X-ray CCD cameras, X-ray Imaging Spectrometer (XIS), which are sensitive from 0.3-12 keV. Data 
were collected from all four XISs as these observations were made while XIS2 was still operational. 
Additionally, there is a Hard X-ray Detector (HXD) which is sensitive from 10-600 keV. It consists of 
Positive Intrinsic Negative (PIN) silicon diodes and Gadolinium Silicate (GSO) phoswitch counters. 
All four observations were made with 3C120 placed at the nominal aim point for the XIS detector. 

Data reduction was performed using the HEASOFT data reduction suite. The event files were 
cleaned using standard filtering criteria and spectra were extracted with the xselect software as 
per the Suzaku ABC guide. The most recent CALDB files at the time of analysis (Fall 2011) were 
used. The 3x3 and 5x5 event files were combined within xselect prior to spectral extraction. The 
cleaned event files contained 144 ks of good data. Spectra were extracted from counts within a 
208" circular region centered on the source. Background spectra were produced from an annular 
region around the source. We generated response files for each XIS sensor and each pointing using 
the xisrmfgen and xissimarfgen scripts. Lastly, for each XIS, the 4 pointings were co-added; 
however we did not co-add spectra from the different XISs. This allows us to remain sensitive to 
cross-calibration issues. 

Similarly, the HXD data was reduced as per the ABC guide. The HXD/PIN was able to 
detect 3C120 while the HXD/GSO was not. We produced cleaned event files and extracted spectra 
using the standard tools. These were then dead time corrected using the hxddtcor script. We 
obtained the appropriate "tuned" non-X-ray background files (NXB) for these observations. Epoch 
appropriate response files were obtained from the Suzaku CALDB website. The Cosmic X-ray 
Background (CXB) files were simulated with XSPEC 12.7 using the provided response file and 
following the typical model based on HEAOl results (Boldt 1987). The CXB and NXB files were 
then added to produce a final background spectra. As with the XIS data, the four HXD/PIN 
observations and their backgrounds were co-added. 
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3. Spectral Analysis 

All of the spectral analysis performed for this paper used XSPEC vl2.7 (Arnaud 1996). Errors 
were computed using Monte Carlo Markov Chain (MCMC) methods, and are reported at the 90% 
level of confidence. The fitting parameters for our spectral model were tied across the HXD/PIN 
and XIS detectors save for a cross-normalization factor to compensate for differences in the absolute 
flux between detectors. The seed value for the HXD/PIN cross-normalization was 1.16 in order 
to be consistent with current recommended calibrations (Suzaku Memo 2008-06) but, as we note 
below, the data strongly requires an HXD/PIN cross-normalization closer to unity 

We begin by examining the spectra in the ranges of 0.7-1.5 keV and 2.5-10 keV. The data in 
the 1.5-2.5 keV band were ignored during all fits due to known issues with the XIS calibration 
over this range of energies. Our initial model consisted of a simple absorbed power-law. Figure 1 
zooms in on the 3-9 keV band and shows that this model results in a poor fit to the data {x 2 /v = 
6556/4768 [1.37]) but illuminates a variety of interesting spectral features. The most noticeable of 
such features is the strong emission line at 6.4 keV in the source rest frame which is likely due to 
the FE-Ka fluorescence line. This line has a red wing which extends to ~5 keV indicating possible 
relativistic broadening effects. There is also a weaker emission feature at ~6.9keV which is likely 
a blend of the cold Fe-K/3 line and the Fe xxvi line. The hard excess flux at energies greater than 
8keV is interpreted as the Compton hump associated with the X-ray reflection producing these 
emission features. 

Guided by this preliminary fit and taking guidance from K07, we built a heuristic model con- 
sisting of photoelectric absorption (phabs), the continuum from reflection off cold matter (pexmon, 
Nandra et al. 2007), a narrow Gaussian Fexxvi-Kq; line at 6.97 keV, and an ionized reflection 
model (ref lionx, Ross & Fabian 2005) convolved with a relativistic smearing model (relconv, 
Dauser et al. 2010) to produce a smeared reflection profile for an ionized accretion disk. The 
pexmon model is a particularly useful evolution of the pexrav model of Magdziarz & Zdziarski 
(1995). In addition to producing the Compton backscattered continuum it also generates physi- 
cally self-consistent Fe-K and Ni-K lines. 

Parameters common to pexmon and ref lionx were tied. Additionally, pexmon takes the low- 
Z metal abundance and iron abundance (Zp e ) relative to solar values as fitting parameters. We 
allowed Zp e to be a free parameter, but the low-Z abundance was frozen at unity. The radial 
dependence of the emissivity of the disk reflection is assumed to be a power-law, e ~ r - ^; broken 
power-law profiles did not lead to a significant increase in the goodness of fit. Previous fits by K07 
and Ogle et al. (2005) have assumed an emissivity index of /3 = 3, whereas we have allowed this 
value to be free. 

We fit this model to the global XIS and HXD/PIN spectra over 0.7-45 keV, considering two 
independent scenarios. Both scenarios are guided by the observation that radio flares are associated 
with X-ray dips (not flares), so we assume X-ray emission is dominated by the disk corona (Ogle et. 
al 2005, Marshall et al. 2009) and not the base of a jet. Scenario #1 involves making the assumption 
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that reflection components terminate at the ISCO and fitting for a. We note the nominal range of 
spin values in relconv is likely not appropriate, with MHD torques from the accretion flow limiting 
it to \a\ < 0.95 (Penna et al. 2010). Thus, we have restricted the spin parameters that we consider 
to the range a G [-0.95,0.95]. This produced a good fit to the data {x 2 /v = 5527/5033 [1.10]). 
We obtained a SMBH spin of a < —0.1, indicating a moderate to rapidly spinning retrograde hole 
(Table 1). We also note a favorable result for the inclination of i = 17.7^g7 degrees. This is in 
good agreement with the constraints derived from super luminal jet motion (i < 14°) and suggests 
that the inner disk and jet are aligned as expected from the theory of Bardeen & Petterson (1975). 

Scenario #2 investigates the possibility that the inner region of the disk is truncated beyond 
the ISCO, an approach similar to the original analysis of K07. With our model this is accomplished 
by freezing the spin at a = 0.95 and allowing the inner disk radius to be a free parameter. The 
choice of a rapid prograde spin permits the model to integrate over the maximum range of allowable 
radii. The fit was nearly statistically identical to the "spin-free" fit {x 2 l v = 5524/5033 [1.10]). The 
fitted inner radii was ri n = ll.7^\r g (Table 1). Essential disk parameters were not preferentially 
seeded with results from the "spin-free" fit; yet they still arrived at very similar values. The full 
set of parameters from both fits are listed in Table 1 (also see Figure 2). 

The strength of the narrow iron line was measured with a separate model consisting of an 
absorbed power law and two gaussian components at 6.4 keV and 6.97 keV. This allows us to isolate 
the narrow Fe-K« line. The equivalent width was found to be EWj<- a = 55 ± 6 eV. Following a 
similar prescription we isolated and measured the strength of the relativistically broadened line 
finding EW Broad = 48 ± 10 eV. We also computed model fluxes and luminosities which can be 
found in Table 2. 

We note that a curious aspect in both of these scenarios is the behavior of the cross-normalization 
factor between the XIS and HXD/PIN datasets. The fiducial value put forth by the Suzaku science 
team, for this dataset, is C = 1.16. Both scenarios achieved a best fit with C ~ 1 which is not in 
statistical agreement with the fiducial value. In order to test the validity of this result we repeated 
our analysis, this time leaving the cross-normalization frozen at C. This produced a statistically 
worse fit to the data, in both scenarios, with Ax 2 ~ 90 for one additional degree of freedom. There 
are strong residuals above 8 keV which indicate an offset normalization for the HXD/PIN. We 
conclude that our model accurately describes the data and the fiducial cross-normalization factor 
is not applicable for this dataset. 

4. Discussion & Conclusion 

We confirm the presence of broad asymmetrical iron line features in 3C120 due to emission 
from the inner regions of an accretion disk. Interpreting the structure of the inner disk from these 
data we arrive at two possible conclusions. Under the assumption that fluorescent iron line emission 
terminates at the ISCO, our analysis of Suzaku data for 3C120 indicates that the central SMBH 
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may be in a state of retrograde spin. The probability distribution for the spin parameter, derived 
from MCMC methods, confirms a spin of a < —0.1 at a 90% confidence. We can reject a rapidly 
rotating prograde black hole with spin a > 0.8 at the 99% confidence level. We also consider 
evidence for a natural truncation of the disk, independent of the spin configuration of the central 
SMBH. In doing so we find that an inner disk radius of r in = ll.T^f^g ^ s strongly preferred by 
our model. Both probability distributions can be seen in Figure 3. 

Theoretically, the possibility of a retrograde spin in 3C120 is intriguing. Employing the flux 
trapping picture of Reynolds, Garofalo, & Begelman (2006), Garofalo (2010) suggested that a 
rapidly-rotating retrograde black hole may be an effective configuration for producing a power- 
ful jet, with the enlarged ISCO of the retrograde hole trapping a large bundle of magnetic flux 
and driving a powerful Blandford-Znajek jet. Interestingly, this idea is not borne out in General 
Relativistic magnetohydrodynamic simulations which find that jets from prograde black holes are 
factors of several more powerful than from the equivalent retrograde hole (Tchekhovskoy et al. 
2012; Tchekhovskoy & McKinney 2012). It is not clear, however, that the "magnetically choked" 
configuration employed by Tchekhovskoy et al. is equivalent to the advective-diffusive equilibrium 
underlying the Garofalo (2010) models and so it is too early to say that such models have been 
invalidated by the simulations. More theoretical work is required to assess the efficiency with which 
retrograde black holes can power jets. 

It is equally intriguing to consider the implications of a natural truncation of the inner disk 
beyond the ISCO. One possibility is the disk undergoes a transition to a radiatively inefficient 
accretion flow (RIAF) at r ~ 10r ff . This would cause the inner regions of the accretion disk 
to become optically thin and extremely hot, truncating iron line emission and all other signs of 
reflection. Such a transition is believed to occur once the accretion rate falls below an Eddington 
ratio (£) of a few percent (Esin et al. 1997). For example, disk truncation is clearly measured in 
MCG -5-23-16 (£ = 0.05; Reeves et al. 2006) and NGC 4593 (£ = 0.04; Markowitz & Reeves 2009). 
However, 3C120 is not a low-£ source. Using our model 2-10 keV luminosity along with a bolometric 
correction of K2-wkeV = 20.6 (Vasudevan & Fabian 2009) and SMBH mass M 3 ci20 = 5.5 x 10 7 M & 
(Peterson et al. 2004) we compute an Eddington ratio of £30120 = 0.23. This is high enough that 
a standard transition to a RIAF at 10r 9 becomes unlikely. 

A more intriguing possibility is that we are seeing a transitory truncation of the disk. Long- 
term monitoring of 3C120 reveals periodic dips in X-ray luminosity coincident with (or slightly 
preceding) the ejection of a new superluminal VLBA knot in the jet accompanied by large radio 
bursts. It has been hypothesized that this may be due to a catastrophic instability, thermal or 
dynamic in nature, which completely destroys/ejects the inner regions of the disk; and that this 
process is an important component of jet formation (Marscher et al. 2002, Chatterjee et al. 2009, 
King et al. 2011). The destruction of the inner line emitting regions of the disk would cause iron 
line emission to cease. 

Is it possible that such a dramatic event in the accretion disk of 3C120 occurred prior to the 
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observations by Suzaku! In fact, a significant X-ray dip followed by a jet ejection event did indeed 
occur at the very beginning of 2006 (2006.00 ±0.03; Chatterjee et al. 2009), just one month prior to 
these Suzaku observations. Following the destruction of the inner accretion disk, the inner regions 
need some amount of time to refill with material; so the issue is whether or not the disk had time 
to reform between the ejection event and the Suzaku observations. At this point, it is instructive to 
consider relevant accretion disk timescales, focusing on a radius of 10 r g while assuming an angular 
momentum transport parameter of a = 0.1 and radiative efficiency of ij = 0.1 for the disk. The 
dynamical timescale is td yn = (r 3 /GM) 1 / 2 ~ 9ks and the viscous timescale is t v i sc = td yn /(a[h/r] 2 ) 
where h/r is the scale height of the disk. Accretion disk theory (Shakura &: Sunyaev 1973) states 
that the inner region of this radiation-pressure dominated accretion disk has a geometric thickness 
of h = (3£/2r?) r g . Using our computed Eddington ratio we find h « 3.45 r g ; which results in 
a viscous timescale of t v i sc ss 750 ks (8.7 days) at 10 r g . Given that the viscous timescale has a 
strong dependence on the precise choice of radius (t v i sc oc r 7 / 2 ), it is reasonable to assume that this 
timescale may be easily extended to a month or more over the range of ri n allowed by our fit. 

We have presented a reanalysis of the 2006 Suzaku data for the AGN 3C120. Thanks to updated 
calibrations and models we were able to reconfirm and expand on prior analysis of this dataset by 
K07. We confirmed the presence of broad iron line features consistent with relativistically-blurred 
ionized reflection from the disk, and show that a disk truncation radius of r ~ 7-14 r g is preferred. 
We arrive at two possible explanations for the truncation; either this black hole is in a state of 
retrograde accretion, or we are seeing the refilling of the accretion disk following a jet-ejection 
event. Detailed spectral monitoring of BLRGs with VLBA in conjunction with XMM-Newton or 
Suzaku are required to assess this last possibility — this should be a key goal of future programs. 

We thank Andy Fabian, Ashley King, Anne Lohfink, Cole Miller, and Margaret Trippe for 
stimulating discussions during the course of this work. We also thank our anonymous referee whose 
comments greatly improved this manuscript. This research has made use of data obtained from 
the High Energy Astrophysics Science Archive Research Center (HEASARC), provided by NASA's 
Goddard Space Flight Center. CSR acknowledges support from NASA through the Suzaku Cycle-5 
GO Program (grant NNX10AR31G) and the ADAP (grant NNX12AE13G). 
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Fig. 1. — The data/model residuals of the simple absorbed power-law fit to the XIS data in the 
3-9 keV range. The iron line complex is clearly visible around 6.4 keV. 
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Fig. 2.— The best-fit model to the global XIS+HXD/PIN spectrum from 0.7-45 keV. The top 
panel shows the individual model components, multiplied by two factors of energy. The bottom 
panel shows the data/model residuals for the "spin- free" fit. The residuals for the "rj n -free fit" are 
nearly identical. 
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Fig. 3. — The probability distributions for a and ri n . The distribution for spin was strongly peaked 
for negative values of a. The choice of plotting log 10 P{a) allows the structure of the distribution to 
remain discernible. The dashed and dot-dashed lines indicate the 90% and 99% CL upper limits, 
respectively. 
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Table 1: Spectral Fitting Parameters 



Model Component 


Parameter 


Suzaku (0.7-45 keV) 


Suzaku (0.7-45 keV) 


Galactic Column 




U.100_ _ooi 


0.1581°;°°} 


PL 


r 


1 07+O.OI 
i -°'-0.01 

(1.241°-°}) x 10- 2 


1-861851 
(i.25i8;8}) x io- 2 


Cold Reflection 


n 


U - OD -0.04 


^7+ - 04 
U -°'-0.05 




PL Cutoff (keV) 


150* 


150* 


Fe xxvi-Ka Line 


E (keV) 


6.97* 


6.97* 




a (keV) 


0.01* 


0.01* 




A^ine 


(7.01^) x 10- 6 


(6.9+}^) x IO" 6 


Accretion Disk 


Zfc 

a 

i 

P 


0.501°;°^ 
201^ 

i 7 7+6.5 
L '-'-8.7 
o qo+0.23 
^• oo -0.30 

ISCO* 


0.501°;°* 
2113 

18.218:2 

o co+0.69 
Z - UO -0.36 

ll.Tiy 




^out 


400* 


400* 


PIN/XIS Normalization 




n qo+0.02 


Q8 + °-° 2 


SMBH Spin 


a 


< -0.1 


0.95* 


X 2 /" 




5527/5033 [1.10] 


5524/5033 [1.10] 



Note — All errors are quoted at the 90% confidence level for one interesting parameter (Ax 2 = 
2.7). Parameters which were fixed for the fit are marked with an asterisk. The units are as 
follows: column density is in 10 22 cm' 2 , normalizations are in photons cm -2 s _1 , iron abundance 
is relative to solar values and linked between pexmon and ref lionx, inclination is in degrees, 
all radii are in r g , and spin is dimensionless (See text for details). 
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Table 2: Model Fluxes and Luminosities for 3C120 



0.7-2 keV 2-10 keV 0.7-10 keV 

Flux^ 1.61 x 10- 11 4.03 x lO" 11 5.64 x lO" 11 

Flux* 2.20 x 10 -11 4.17 x 10 -11 6.39 x 10 -11 

Luminosity 3.95 x 10 43 1.00 x 10 44 1.40 x 10 44 

Luminosity^ 5.52 x 10 43 1.05 x 10 44 1.60 x 10 44 



Note — The model fluxes and luminosities for 3C120. Fluxes are given in erg cm~ 2 s -1 . 
Luminosities are computed with a standard cosmological model of Hq = 70 km s _1 Mpc -1 , 
VL\ = 0.73 and Qm = 0.28 and are given in erg s _1 . Values marked with a f have been 
absorption corrected. The errors are dominated by systematic effects and thus not reported. 



